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Abstract 

In condensed matter physics, nonequilibrium pump-probe time domain spectroscopies have be- 
come an important tool to disentangle degrees of freedom whose coupling leads to broad structures 
in the frequency domain. Here, using the time-resolved solution of a model photoexcited electron- 
phonon system we show that the relaxational dynamics are governed by the equilibrium self-energy 
so that the phonon frequency sets a window for "slow" versus "fast" recovery. The overall tempo- 
ral structure of this relaxation spectroscopy allows for a reliable and quantitative extraction of the 
electron-phonon coupling strength without resorting to an effective temperature model. 



In pump-probe experiments, phonons provide the primary channel for photoexcited elec- 
trons to relax and reduce their energy via scattering. By emitting phonons, the electrons 
dissipate the excess energy imparted by the photons, eventually reaching a new thermal 
equilibrium. The rate at which electrons can dissipate their energy is determined, in part, 
by the strength of electron-phonon coupling. This coupling is also of critical importance in 
determining other properties of a material, such as electrical or thermal transport, the effec- 
tive electronic mass enhancement, the size of the polaronic cloud that surrounds an electron 
in the normal state, or the superconducting properties of a BCS superconductor. As with 
many basic coupling constants in nature, it is difficult to directly measure the strength of 
electron-phonon coupling, making it a challenge to assess its role in materials like high tem- 
perature superconductors. Here we show how one can measure the electron-phonon coupling 
from relaxational pump-probe photoemission spectroscopy. 

The strength of electron-phonon coupling is commonly described through one of two 
dimensionless coupling constants: A determines the change in the electronic specific heat 
due to electron-phonon coupling, the renormalization of the Fermi liquid properties due 
to phonons, and the superconducting properties, while the transport coupling strength X tr 
determines transport properties. In conventional superconductors, an inversion 1 of the fine 
structure in the tunneling density of states above the gap can be used to determine A to high 
precision. The transport coupling strength \ tr can be found by inverting optical conductivity 
data, 2,3 even in the normal state. 4 In both cases, the phonon effects must be dominant and 
strong for the inversion to be accurate. Other experimental methods of determining electron- 
phonon couplings include x-ray scattering, 5 Raman scattering, 6 or neutron scattering. 7 

Angle-resolved photoemission (ARPES) has been widely used to determine the low-energy 
electronic structure of a variety of materials. 8,9 Many-body interaction effects are manifest 
via the self energy E giving rise to renormalizations of the electron mass and an overall 
broadening of the spectral function. As an example, it was thought that ARPES data for 
the cuprates could thereby be studied to extract the boson coupling that is responsible 
for high temperature superconductivity. 10-26 These investigations have provided important 
insights into the underlying bosons and interactions responsible for kinks in the electronic 
dispersion. However, extracting the self-energy and a meaningful coupling constant from 
frequency-dependent ARPES data can be challenging, especially when many sources of 
incoherence are present simultaneously. 
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As an alternative, a time domain approach using the relaxation of excited electrons back 
to equilibrium as a direct technique to measure electron-phonon coupling was proposed by 
Allen. 27 The basic quantities for this method are the deviations from the initial equilibrium 
values of the electron and phonon occupations in momentum space. A semi-classical Boltz- 
mann approximation that averages over momenta was used to extract estimates for electron- 
phonon couplings from pump-probe photoemission data in simple elemental metals. 28 ' 29 It 
also was applied to high-temperature superconductors 30-33 even though the approximation 
can be justified at best for simple Drude metals subject to relatively weak driving fields in 
the linear response regime. 

However, in more correlated materials, such as transition metal oxides like the cuprates, 
or in multiband systems with multiple disconnected Fermi surface pieces, such as the iron 
pnictides, the momentum dependence must be resolved, since it is a crucial ingredient for 
these materials' interesting properties. 34,35 In addition, a full quantum mechanical formula- 
tion is required to treat arbitrarily strong fields and strong energy and momentum dependent 
scattering, particularly in materials where the balance between various competing phases 
may be tipped with the application of strong electric fields. Examples include possible 
light-induced superconductivity in a stripe-ordered cuprate, 36 a dramatically enhanced con- 
ductivity in an oxide heterostructure. 37 , coherent phonon oscillations serving as a marker 
for a photoinduced phase transition in VO2, 38 , and pump-probe photoemission providing a 
novel classification scheme for charge-density wave insulators in the time domain. 39 

In this paper we present a microscopic quantum nonequilibrium approach with full mo- 
mentum and energy resolution that can be applied when the semiclassical Boltzmann equa- 
tion fails, and can be used in conjunction with standard many-body and material specific 
computational approaches. Using our model system, we show how the time (t), momentum 
(fc), and energy (u) resolution provided by state-of-the-art tr-ARPES 30-33 ' 40-45 experiments 
in conjunction with the dependence of extracted decay times on the initial equilibrium sam- 
ple temperature provides a direct method to measure the dominant phonon energy and 
electron-phonon coupling. Hence, we show that relaxational spectroscopy in the time do- 
main allows the possibility to extract detailed information about the system in equilibrium 
and lend itself to a simple interpretation beyond the traditional two-temperature model 
phenomenology. 
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I. RESULTS 



A. Photoexcited electronic structure: Energy and temperature dependent relax- 
ation 

We consider 2D tight binding electrons in a single partially filled metallic band cou- 
pled to the lattice via electron-phonon coupling in a Holstein model 46 (see Supplementary 
Material). Phonons are described by a single optical branch with momentum-independent 
phonon frequency f2 = 0.1 eV and electron-phonon coupling g = ^0.02 eV as representative 
values for optical phonons in transition metal oxides. Our approach can be generalized to 
momentum-dependent phonons and momentum-dependent couplings. In equilibrium, the 
phonon frequency sets the relevant energy scale at which a kink at ±fi in the band structure 
occurs 47 , as shown in Figure la-c. The dimensionless coupling constant A is a measure of 
the strength of lattice distortion relative to the kinetic energy of the electrons. A is defined 
via the slope of the real part of the self-energy and determines the electronic mass renor- 
malization. Our choice of parameters corresponds to a value of A = 0.42. This value of A 
is sufficiently large to produce a kink in the dispersion, but still small enough to justify the 
use of the Migdal approximation for the electron-phonon coupling (see Supplement). 

The system is prepared in an initial equilibrium state at a temperature T. Subsequently 
the electrons are photoexcited by a pump laser pulse modeled by a time-dependent electric 
field. We utilize the Keldysh technique for calculating the electron propagation along the 
Keldysh contour (see Supplement for details). The electric field of the pump pulse is included 
via a Peierls substitution valid for any temporal field shape or strength. This approach has 
been formulated earlier by Jauho and Wilkins 48 which we extend numerically to longer time 
scales and arbitrary field strengths. Specifically we choose a Gaussian field envelope of width 
a = 25 femtoseconds (fs) and a frequency of 0.5 eV, as representative of a typical pump 
setup. The delay time is defined with respect to the center of the Gaussian pulse envelope. 
The electric field is linearly polarized along the Brillouin zone diagonal (see inset in Figure 
la and Supplement). We have checked in the analysis of the transient relaxation that details 
of the pump are not especially important unless there is a competing phase to which we can 
pump, 36 ' 49 and unless we are interested in the behavior during the pump. The subsequent 
tr-ARPES probe pulse has a Gaussian envelope width of a w = 103 fs. In the Supplement, 
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we also show the equilibrium tr-ARPES spectra for smaller a pr , which enables a better time 
resolution but necessarily leads to worse energy resolution. 

In Figure 1 we show the energy- and momentum-resolved tr-ARPES spectra along a nodal 
cut for three delay times before (la), during (lb), and after the pump pulse (lc). The pump 
pulse drives electrons from below the Fermi level Ep to empty states above Ep (Figure lb) as 
the system absorbs energy leaving the electrons in a photoexcited state. The excess energy 
of the electrons is then transferred into the bath of phonons, and the electrons relax back 
towards equilibrium (Figure lc). Figure Id provides a comparison of the fe-resolved energy 
distribution curves (EDC) before and 558 fs after the pump. Clearly, the phonon window W 
= [-Q, Q] around the Fermi level sets an energy scale for persistent spectral changes: while 
electronic quasiparticles outside this window have already relaxed back to their equilibrium 
EDC, the quasiparticles inside remain excited at the particular time delay shown here. This 
phonon window effect is highlighted in the momentum-integrated pump-induced spectral 
changes in Figure le: photoexcited electrons above Ep and holes below Ep only persist 
inside the window set by the phonon frequency Q. 

In equilibrium at zero temperature, kinematics prevents electron relaxation unless the 
quasiparticles have sufficient energy to emit phonons, which is why this window is not 
unexpected. In Figure If, we show that raising the initial equilibrium temperature leads to 
smaller spectral changes compared to Figure Id at the same time delay. Correspondingly, 
the fe-integrated spectral changes shown in Figure lg are significantly smaller than the ones 
in le (shown on the same vertical scale), although the phonon window effect is still visible. 

A simple way of thinking about pump-probe experiments has been that pumping the 
system is approximately the same as heating. Models of "hot electrons" assume that the 
laser pump initially raises the electronic temperature. 27 ' 50 The electrons then transfer their 
excess energy to the lattice until the electronic and lattice subsystems equilibrate. Here, the 
hot electron picture does not provide a valid description of the photoexcited electrons in our 
system, which is demonstrated in Figure 2. 

In Figure 2a and 2b, we show the momentum-resolved tr-ARPES spectrum along the 
diagonal momentum cut at a short time delay (4 fs) for two different pump pulses. Figure 
2a is for the same maximum field strength i? max = 0.4 V/a that was used in Figure 1; Figure 
2b is for a ten times stronger field. The weaker pump mainly excites electrons within the 
bare energy band, but the stronger pump clearly disrupts the band structure and induces 
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ladder-like structures which are reminiscent of Wannier-Stark ladders in strong dc electric 
fields. 51 ' 52 

The breakdown of the "hot electron" picture becomes more evident when we integrate 
the tr-ARPES spectra in the momentum window around the Fermi momentum between k\ 
and (see Figure 1). The energy distribution curves in Figure 2c are not simply featureless 
decaying exponentials, as they would be if the "hot electron" model was valid, but show 
distinct higher energy peaks related to the pump characteristics. Although it is possible 
to extract an effective temperature from the linear slope of the logarithmic intensities in a 
very small energy window (< 0.1 eV) near the Fermi level (see Figure S3 in Supplement), 
this "temperature" does not provide a meaningful description of the photoexcited system. 
In particular, it fails to account for the actual excess energy in the electronic subsystem 
due to the spectral features at higher energies. These nonthermal features are even more 
pronounced for the strong pump. This observation clearly shows that photoexcitation is not 
the same as heating, and that a full nonequilibrium modeling beyond effective temperature 
phenomenology is essential for correctly describing pump-probe spectroscopies. 

B. Extracting equilibrium system properties by pump-probe photoemission 

As an effective temperature model is inadequate for our purposes, we now focus on a 
quantitative analysis of energy-dependent relaxation rates for the weaker pump and the 
higher initial temperature T = 290 K. Figure 3a shows the pump-induced spectral changes 
relative to the initial state, integrated along momenta on the Brillouin zone diagonal, in a 
continuous map as a function of energy and delay time. Here, the phonon window effect 
becomes even more obvious: persistent excitations of electrons and holes are only observed 
inside W, whereas photoexcited electrons outside W decay on a much faster time scale. The 
dashed lines mark the phonon window W, which clearly sets the relevant energy scale for 
slow versus fast relaxation. In Figure 3b, we present cuts for selected energies of the same 
data as in Figure 3a: Both electrons and holes inside W decay on a significantly longer time 
scale than their counterparts outside W. The slight asymmetry between holes and electrons 
(negative and positive binding energies) is related to asymmetry in the underlying band 
structure, made clear in the discussion which follows. 

We fit decaying exponentials (solid lines in Figure 3b) to these curves starting at a delay 
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time of 250 fs, chosen to be sufficiently large such that pump and probe pulses do not 
significantly overlap. From these fits we extract the tr-ARPES relaxation rates shown in 
Figure 4a on a logarithmic scale. The relaxation rates show a steep decline inside the phonon 
window. This characteristic behavior indeed is reminiscent of the equilibrium scattering rates 
(self-energy) due to coupling between electrons and an optical phonon branch, reflected in a 
suppression of the imaginary part of the equilibrium self-energy in the phonon window. To 
make this connection between tr-ARPES relaxation rates and the equilibrium self-energy 
apparent, we plot —2 Im S aver in Figure 4a. For comparison we account for energy resolution 
due to the probe pulse width by an additional averaging of relaxation times appropriate for 
the pump-probe setup (see Supplement). 

The good agreement between the tr-ARPES relaxation rates and the equilibrium self- 
energy suggests that relaxational pump-probe photoemission spectroscopy can be used as an 
effective tool for accurately measuring equilibrium system properties. Traditionally, equilib- 
rium ARPES experiments are used to extract an estimate for electronic lifetimes through the 
linewidths of momentum distribution curves (MDCs). For optimal energy and momentum 
resolution and only a single source of decoherence, the full width at half maximum (FWHM) 
is related to the relaxation rate at the corresponding energy via t~ 1 {uj) = FWHM(u)vF/h, 
where vp is the bare Fermi velocity along the corresponding momentum cut and the equation 
is valid near the Fermi level for a linear band dispersion. In order to compare relaxation 
rates extracted from MDC widths directly for the same energy resolution as the tr-ARPES 
decay rates, we extract the MDC widths from a tr-ARPES snapshot using the same probe 
width as for the pump-probe simulation. 

Even in the absence of other line broadening channels such as impurity or electron- 
electron scattering, the restricted energy resolution due to the temporal width of the probe 
pulse leads to MDC linewidths which are not in accordance with the corresponding tr- 
ARPES relaxation rates (see Figure 4a). The MDCs linewidths imply larger rates, which 
becomes apparent inside W where the relaxation rates vary strongly with binding energy 
and the energy resolution therefore has a crucial effect. This demonstrates that pump-probe 
ARPES in principle allows one to extract the equilibrium electron-phonon scattering rates 
more accurately than ARPES MDC measurements. 

As a consequence, more detailed knowledge about the electron-phonon scattering pro- 
cesses allows us to quantitatively extract the electron-phonon mode coupling strength from 
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a comparison of the tr-ARPES relaxation rates T tr (co) and the self-energy E(o>). To this 
end, we treat our tr-ARPES relaxation rates as "experimental data" and assume that we 
know the line shape of Im E(o;). Within the Migdal approximation to the Holstein model, 
the relaxation rate is proportional to g 2 and of the form T(cj) = g 2 4>{ui) (see Supplement). 
The function <p{u) reflects electron-phonon scattering processes and takes into account the 
available scattering phase space related to the electronic band structure via the local elec- 
tronic density of states N{oj). This is the reason why an asymmetry in the band structure 
results in an electron-hole asymmetry in the relaxation rates. 

The measured tr-ARPES relaxation rates r tr are naturally convolved with the probe pulse 
shape and therefore include an additional energy averaging, since there is a tradeoff between 
time resolution and energy resolution. Therefore we apply an energy averaging procedure 
to 4>{uj) (see Supplement), which yields ave r- The electron-phonon coupling strength is then 
extracted from g 2 = r tr (a>) / '0 aver (co'). Figure 4b shows this ratio as a function of energy, 
which serves as an a posteriori estimate of g 2 given the numerical simulation data. The 
extracted mean value including an error estimate from the statistical standard deviation is 
g 2 = 0.021 ± 0.004 eV 2 , to be compared with the exact input value g 2 = 0.02 eV 2 . This 
agreement highlights the fact that a nonequilibrium measurement can serve as an accurate 
tool for extracting equilibrium system properties, a feat which may be difficult to accomplish 
by pure equilibrium measurements. 

It might seem counterintuitive that a system out of equilibrium relaxes at the equilib- 
rium scattering time scale, determined by the imaginary part of the self-energy, but this 
result should hold in nearly all systems at long times. This occurs as for small deviations 
from equilibrium the change of the spectrum is second order in the deviation from equilib- 
rium, while the change of the occupation is first order, primarily due to the change of the 
distribution function from the equilibrium Fermi-Dirac distribution. 53,54 Hence, the Green 
functions should be the same as in equilibrium, except for a small change of the distribution 
function. If the self-energy of the system in equilibrium does not depend too strongly on 
temperature, then it is weakly dependent on the distribution function, and as a result, we 
will see the nonequilibrium system relax back to equilibrium at the longer time scales via 
the equilibrium relaxation rates. 
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II. DISCUSSION 



The central result of this paper is that equilibrium ideas can be used to understand the 
nonequilibrium dynamics of a pumped system. This is facilitated by the connection between 
relaxation rates for the decay of photoexcited electrons and the equilibrium electron-phonon 
scattering rates. Importantly, both the dependence of the rates on the binding energy and 
on the initial equilibrium temperature affect the result, but not the conclusion. We have 
shown that a simplified modeling in terms of effective temperatures does not provide an 
adequate description of photoexcited systems with full nonequilibrium quantum dynamics. 
However, the relation between the equilibrium self-energy and the nonequilibrium relaxation 
makes pump-probe spectroscopies a convenient tool to disentangle equilibrium decoherence 
pathways. As a consequence, the strength of electron-phonon couplings can be extracted 
from pump-probe relaxation times with a high level of accuracy. In addition, our approach 
also can be applied to general electron-boson interactions, for instance coupling to magnons. 
It could therefore potentially be used to better characterize the bosonic couplings in high- 
temperature superconductors. 



III. METHODS 

We solve a set of integro-differential equations 55 along the time contour in order to ac- 
count for electron-phonon interaction effects via a double time self-energy, which dresses the 
noninteracting Green function. The electric field of the pump pulse is included via a standard 
Peierls substitution. The code is parallelized easily since the equations for different momen- 
tum points decouple in the discretized Brillouin zone. An additional Fourier transform of 
the lesser component of the resulting interacting Green function is performed to extract the 
tr-ARPES spectra. 56 In this step, a shape function is introduced to take the probe pulse 
width into account. Additionally, the momenta need to be time-shifted depending on the 
pump field to ensure gauge invariance. 57 The details of the formalism are presented in the 
Supplement. 
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FIG. 1. fc-resolved photoexcited electronic structure, (a) Equilibrium tr-ARPES intensity 
at temperature T = 116 K for a cut along the Brillouin zone diagonal. The pump field polarization 
is aligned along this direction (inset). Horizontal dashed lines mark the position of "kinks" in the 
electronic dispersion at ±Q. (b) Map of the tr-ARPES intensity shortly (87 fs) after the pump 
pulse and (c) at a later delay time t = 558 fs. The system has partially relaxed toward equilibrium 
outside the energy range ±0, but quasiparticles remain excited within this range, (d) Energy 
distribution curves (EDCs) corresponding to the data shown in (a, equilibrium) and (c, 558 fs) for 
momenta k\ to as indicated in panels (a) and (c). For clarity EDCs are vertically shifted, (e) 
Intensity difference AI U between photoexcited (558 fs) and equilibrium tr-ARPES data. (f),(g) 
EDCs and spectral changes (same as in (d),(e)) for a higher initial equilibrium temperature T = 
290 K, showing that the time for recovery shortens for increased temperature. 
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FIG. 2. Weak versus strong pump field, a, tr-ARPES intensity map (delay time of 4 fs) 
for a weak pump field (same as in Figure 1) with maximum field strength £ max = 0.4 V/arj. b, 
Ten times stronger pump field (E m8iX = 4 V/ao). For a lattice constant ao = 3.4 A, these peak 
field strengths are ~ 12 MV/cm ("weak pump") and 118 MV/cm ("strong pump"), respectively, 
c, Energy distribution curves integrated over momenta along the Brillouin zone diagonal (same 
momentum window k\ to % close to kp as indicated in Figure 1) for various pump-probe delay 
times. 
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FIG. 3. Transient relaxation of the pump-induced spectral intensity changes for equi- 
librium temperature T = 290 K. (a) Spectral changes as a function of pump-probe delay time 
integrated over momenta along the whole Brillouin zone diagonal, (b) Same data as in a for se- 
lected energies inside (0.05 eV) and outside (0.25 eV) the phonon window for positive (electrons) 
and negative (holes) binding energies. 
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FIG. 4. Energy dependent relaxation due to electron-phonon scattering, (a) Relaxation 
rates from tr-ARPES (red circles) compared with equilibrium MDC widths (green diamonds) for 
the same probe pulse width. The orange line shows the corresponding equilibrium electron-phonon 
scattering rate taking into account energy resolution (see Supplement), (b) Extraction of the 
electron-phonon coupling strength from the ratio of tr-ARPES decay rates and equilibrium self- 
energy (see Supplement). The red solid line shows the mean value of the data points. The black 
dashed line indicates the exact input value of g 2 . The grey shaded area shows the confidence 
interval estimated from the statistical standard deviation of the data. 
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SI. FORMALISM 



For the formulation of the nonequilibrium problem, we employ the well-established 
Keldysh Green function technique including the electron-phonon coupling via the electronic 
self-energy in the Migdal approximation. 1 The double-time Green functions allow us to ex- 
tract the relevant information (occupation of states, tr-ARPES response) about the excited 
electron system in a microscopic way. The real and imaginary times involved in the problem 
define the Keldysh contour shown in Figure SI. The resulting integro-differential equations 
of motion for the Green functions are solved using massively parallel numerical integration. 

We study the Holstein Hamiltonian 

H = £ e{k)c{c k + n b l b i " 9 E + &!), (SI) 

k i i 

where the individual terms are the electronic kinetic energy with dispersion e(k), the phonon 
energy for dispersionless Einstein phonons of frequency Q, and the electron-phonon interac- 
tion, respectively. Here is the electronic annihilation (creation) operator in momentum 
space (real space), bf^ is the bosonic phonon annihilation (creation) operator on lattice site 
i, and g is the momentum-independent electron-phonon coupling strength. 

The driving field is included directly in the propagator via Peierls substitution in the 
double-time formalism, which leads to the bare non-equilibrium Green function 2 



G° k (t,t') =i [n F (e(k))-e c (t,t')\ 



x exp 



-if dt e{k- A(t)) 
Jt> 



(S2) 



where n F (u) = 1/(1 + exp(w/T)) is the Fermi function at temperature T, t and t' are times 
on the Keldysh contour C, 9 c (t, t') — 1 if t is later on the contour than t' and otherwise, and 
A(t) is the time-dependent vector potential of the driving field. e(k) is the single-particle 
energy dispersion for a square lattice with nearest (V nn ) and next nearest neighbor (V nnn ) 
hopping, 

e(k) = —2V nn (cos k x + cos k y ) + AV nnn cos k x cos k y — pi, (S3) 

measured with respect to the chemical potential /i. In the following we will use a fixed set 
of model parameters V nn = 0.25 eV, V nnn = 0.075 eV, \x = —0.255 eV. A(t) is the vector 
potential at time t, which is related to the electric field via A(t) = — f f E(t')dt'. Here, we 



use the convention that Ti — c — e — 1, and we work in the Hamiltonian gauge, i.e. the 
scalar potential is set to zero. Energies are measured in units of eV, time scales in units of 
h/eV « 4.1 fs. 

We incorporate the electron-phonon interactions in the Migdal limit, which is appropriate 
for weak coupling. In this case the electronic self-energy due to electron-phonon coupling is 

S(t,0 = <^E £, °(*'OG*(*,0- (S4) 

The bare phonon Green's function D°(t,t') is 3 

D°(t, t')=-i [n B (ty + 1 - 6 c (t, t')} exp (iQ(t - t')) 

- i [n B (H) + 9 c (t, t')] exp (-iQ(t - t')) (S5) 

where = l/(exp(f2/T) — 1) is the Bose function. 

With the self-energy above, we solve the Dyson equation for the interacting Green func- 
tion G, 

G k (t, t') = G° k (t, t') - j dt l( lt 2 G k (t, iOS^i, h)G k (t 2 , t'). (S6) 

This can be done by casting the Dyson equation as a matrix equation. As an alternative, 
the Dyson equation can be rewritten by expanding the integral with Langreth rules, which 
gives equations of motion for the Matsubara, retarded, real-imaginary, and lesser Green's 
functions. 4,5 This leads to a set of Volterra integro-differential equations that can be solved 
via standard numerical integration. For the results presented here, we have used the Volterra 
scheme. Its main advantages are that the equations are manifestly causal, and the required 
memory is reduced by the use of symmetries. 

The external electric pump pulse field is included via a standard Peierls substitution. 6,7 
The pulse that is of direct interest to pump-probe experiments is, by nature, a propagating 
light pulse; this implies an oscillating field without a zero-frequency component. We model 
the pump via an oscillating vector potential along the (11) or Brillouin zone diagonal direc- 
tion with a Gaussian profile, where we assume that the field is slowly varying spatially and 
thus neglect the spatial dependence: 

(t — t ) 2 

A(t) = (x + £)F max sin(cV) exp ( - 1 ° ; ) . (S7) 
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The tr-ARPES intensity for a probe pulse of width a pr is computed from 8 



I{k, uj, t ) =Im 



x e 



-(t-t a fl2al Te -(t'-U,Yl2al Te ^{t-t') 



(S8) 



where the integral is along real times, and each component fcj of the momentum argument of 



of the vector potential) to restore gauge invariance. 2 

Figure S2 shows the probe pulse width dependence of tr-ARPES intensity maps and 
corresponding EDCs, revealing that the necessary energy resolution to resolve the phonon 
kink and the sharpening of spectral changes to be inside the phonon range has to be achieved 
by trading off intrinsic time resolution. 



S2. BREAKDOWN OF EFFECTIVE TEMPERATURE DESCRIPTION 

Figure S3 shows fits of Fermi functions with effective temperatures and chemical poten- 
tials to the low-energy structure of the tr-ARPES energy distribution curves. This figure 
serves to show that one can make such a fit in a small energy window and extract decay 
constants for the time-dependent temperatures and chemical potentials. However, these do 
not provide an adequate description of the excited system, which is not at thermal equi- 
librium. Notably, the time-dependent spectra deviate considerably from the fits already at 
slightly higher binding energies than the fit window. The figure therefore demonstrates the 
breakdown of an effective temperature description and calls for a more advanced data anal- 
ysis. Additionally, we note that even if a Fermi function fit seems to work (see, e.g., Ref. 9), 
it is easily shown that the extracted decay time for the effective temperature is not related 
in a simple fashion to energy-dependent decay times for the spectral intensities. This is due 
to the fact that the intensity at a given binding energy depends nonlinearly on the effective 
temperature in this case, namely through a Fermi function. Therefore it is doubtful which 
information can be gained from the temperature decay time. 



G%(t,t') is shifted (fc; fc; 



1 

t-t> 



J* dtAi(t), where Ai is the corresponding i-th component 
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S3. RELAXATION RATES 



In Figure 1 of the main text, we showed how the relaxation process inside the phonon 
window is strongly temperature dependent. For low temperatures, the time scales become 
considerably longer and thus lead to a slower relaxation. Figure S4 shows the correspond- 
ing equilibrium relaxation times. Obviously, the times due to electron-phonon scattering 
are basically independent of temperature outside the phonon window. Inside the phonon 
window, however, they are strongly temperature dependent. The time scales for our model 
parameters are well below a picosecond for the higher equilibrium temperature T =290 K, 
but become as large as tens of picoseconds for T = 116 K. 

The relaxation times shown in Figure S4 are obtained from the imaginary part of the 
equilibrium retarded self-energy according to 

^ - 2T^W < S9 » 

where 

Im E fl (w) = -ng 2 [N(w + Q) (n B (Q) + n F (n + u)) 

+N(u -n)(n B (H) +n F (Q-u))], (S10) 

and N(ui) = L^ 1 J2 k 5{u — e(k)) (the sum is over the two-dimensional Brillouin zone, with L 
the number of momentum points summed over) is the bare electronic density of states. 10 The 
relaxation rate T(ui) = l/r(w) is thus of the form g 2 <p(u), which is used for the extraction 
of g 2 from the tr-ARPES data in the main text. Here, the equilibrium self-energy is taken 
at the initial temperature of the pump-probe simulation. 

For the relaxation rates shown in Figure 4 of the main text, we apply an additional 
averaging procedure to account for the energy resolution due to the finite width of the probe 
pulse. The averaged self-energy is obtained by convoluting its inverse with a resolution 
function: 

1 f 1 f —uj' 2 \ 1 

Im S avcr (u;) = J dUJ '^V^ 6XP {^iJ ImEV)' (S11) 

For the probe pulse width a pr = 103 fs, the corresponding energy resolution is given by a TCS 
= V2h/a pr w 57 meV. Note that this is in contrast to Matthiessen's rule for scattering times, 
which states that inverses of times should be averaged instead of the times themselves. The 



difference is that Matthiessen's rule applies for averages of scattering processes for various 
sources of scattering, all for the same particle. Here, however, we are interested in averages 
of time scales for different particles. The reason why one should average the relaxation 
times here, instead of the rates, is the following: In the long-time limit, the slowest decaying 
components determine the effective relaxation rate. To give those components with the 
longest time scales the largest weight, one has to average times, not rates. The same 
averaging procedure is applied to yield (f) aver , which is used to extract g 2 in Figure 4b of 
the main text. 

Finally, we demonstrate that the relaxation times for particular energies do not signifi- 
cantly depend on the pump field strength as long as one does not measure times too close 
to the application of the pump. To this end, we show in Figure S5 a comparison of the time 
scales for exponential fits to the time traces of the energy-dependent tr-ARPES intensities 
for the weak and strong pumps (see Figure 2 of main text). The times extracted for weak 
and strong pump pulses agree reasonably well. This overall agreement of the time scales 
demonstrates that the interpretation of the nonequilibrium dynamics in terms of equilibrium 
electron-phonon scattering rates is valid even beyond the weak pump regime. 
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FIG. SI. Pump-probe photoemission setup with Keldysh contour. The time-resolved 
pump-probe photoemission is modeled using nonequilibrium Green functions defined on the 
Keldysh contour. The contour runs from the initial time i m j n to a maximum time £max cover- 
ing the time span of the pump pulse (red) of temporal width a and the probe pulse (blue, only 
Gaussian envelope shown here) of width a pr , then back to t m i n , and from there along the imaginary 
time direction to i m j n — i/T. The tr-ARPES snapshots show the electronic structure in equilibrium, 
in the transient regime at short delay times, and in the long-time limit. 
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FIG. S2. Sharpening around phonon edges for increasing energy resolution. The tr- 

ARPES signals for various probe pulse widths become sharper for broader probe pulses. At the 
same time, the restriction of slowly decaying spectral regions to the energy range inside the optical 
phonon range becomes more apparent. 
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FIG. S3. Effective temperatures, (a) Momentum-integrated tr-ARPES data for weak pump 
(see Figure 2 of main text) with fits to Fermi functions A/(ex.p((uj — /x e fl(i))/(^B^eff(*))) + 1) f° r 
time-dependent effective temperatures T e g(t) and chemical potentials fi e s(t). in the energy window 
between and 0.1 eV. (b) Extracted fJ> e s(t) and T e g(t) as a function of delay time. The solid lines 
are fits to decaying exponentials starting at a delay time of 80 fs with decay constants as indicated. 
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FIG. S4. Equilibrium temperature dependent relaxation times. The relaxation times given 
by the self-energy (without additional frequency broadening) for the two equilibrium temperatures 
(116 K and 290 K) used in this paper. 
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FIG. S5. Comparison of tr-ARPES relaxation times for weak and strong pumps. For 

the weak and strong pump fields (same as in Figure 2 of main text), the energy dependence of 
tr-ARPES relaxation times is shown here. 
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